Role of electric field in surface electron dynamics above the vacuum level 
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Scanning tunneling spectroscopy (STS) is used to study the dynamics of hot electrons trapped 
on a Cu(100) surface in field emission resonances (FER) above the vacuum level. Differential 
conductance maps show isotropic electron interference wave patterns around defects whenever their 
energy lies within a surface projected band gap. Their Fourier analysis reveals a broad wave vector 
distribution, interpreted as due to the lateral acceleration of hot electrons in the inhomogeneous tip- 
induced potential. A line-shape analysis of the characteristic constant-current conductance spectra 
permits to establish the relation between apparent width of peaks and intrinsic line-width of FERs, 
as well as the identification of the different broadening mechanisms. 

PACS numbers: 73.20.At, 68.37.Ef, 71.20.Be 
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INTRODUCTION 

A detailed knowledge of the electron dynamics at sur- 
faces is crucial for an understanding of a large variety 
of processes, ranging from electron scattering at surfaces 
to charge transport dynamics across interfaces, relevant 
to design electronic devices 0, 0] . Electrons trapped in 
unoccupied long-lived resonances represent an interest- 
ing workbench. They favour the localization of photo- 
injected electrons at molecular resonances, thus enhanc- 
ing the catalytic activity of metals Q • They also repre- 
sent a valuable probe to investigate the rich phenomenol- 
ogy behind charge injection and hot-electron quenching. 
Experimental techniques such as inverse photoemission 
two-photon photoemission [f|, or ballistic electron 
scattering Q have been traditionally used to study hot 
electron dynamics at surfaces. 

Recently, scanning tunneling spectroscopy (STS) has 
proved to be a useful tool to provide quantitative infor- 
mation about the electronic structure [j], la [9|, |l0|, |l2j and 
also the electron and hole dynamics at metal surfaces 



scale imaging of diamond [25| . A promising application 
of FERs is to provide information about the dynamics of 
electrons in image states at surfaces 15( • This is intrigu- 



II El El- In most cases, these studies have been 



restricted to low applied bias voltages, where the applied 
electric field does not play an important role [181 ]. A 
renewed interest has emerged in using the STM in the 
field-emission regime, i.e., at bias voltages larger than 
the tip work function. In this regime, the applied elec- 
tric field lifts up the potential barrier above the vacuum 
level of the sample, introducing a new class of resonances 
that are absent at low bias voltages, the so-called field- 
emission resonances (FER) [lj| [5(|. In previous stud- 
ies, FERs have been used to explore local changes of the 
surface work function 12 ll 1 221 . scattering properties of 
surfaces and interfaces 123, \24i, and to achieve atomic- 



ing since FERs are a characteristic of the tip-induced po- 
tential barrier itself and, therefore, they would exist even 
in the absence of an image potential. Hence, a model is 
needed which describes the dependence of field emitted 
electrons dynamics along the surface on the topology of 
the surface potential and accounts for STS spectra in a 
wide sample bias range. 



In this paper we demonstrate that electrons trapped 
in long-lived FERs are sensitive to the potential gradient 
induced by the STM tip along the surface. Scattering of 
quasi-free FERs electrons with surface defects gives rise 
to isotropic two-dimensional (2-D) wave patterns, whose 
wave vector components in reciprocal space reflect the 
local perturbation of the surface image potential induced 
by the STM tip. A combined theoretical and experimen- 
tal analysis of the FERs peaks in dl/dV spectra reveals 
that their line-shape carries information about the scat- 
tering properties of the surface, and hence, about their 
band structure. Our calculations permit to identify the 
different intrinsic and extrinsic broadening mechanisms 
of peaks associated with FERs in conductance spectra in 
a wide energy range. The organization of the paper is as 
follows. Section describes the way the experiments were 
done. In section we present the results and discussion 
of them in two subsections: subsection is devoted to the 
analysis of wave patterns that appear in dl/dV maps, 
while subsection presents a line-shape analysis of dl/dV 
spectra. Finally, in section the conclusions of our work 
are presented. 
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EXPERIMENT 

The experiments were performed in a custom made 
ultra-high vacuum scanning tunneling micro scop e in 
thermal equilibrium with a liquid helium bath [261 ] . All 
spectroscopy data presented in this work were acquired 
at 4 K. The Cu(100) sample surfaces were cleaned by 
repetitive cycles of Ar + sputtering (IkeV) and anneal- 
ing at 900 K. The differential conductance (dl/dV) was 
measured using a lock-in amplifier above the low-pass fre- 
quency of the feed back loop (f QC ~3 kHz). The dl/dV 
spectra shown here are taken in constant current mode 
(feed back loop closed) [ljj [2(| • 



RESULTS AND DISCUSION 
Wave patterns 

The electron dynamics of FER states is essentially 
quasi- free in the plane parallel to the sample surface, 
since the corresponding wave functions lie mainly on the 
vacuum side of the surface. Therefore, these electron 
states, similar to image-states, are not affected by the 
corrugation of the surface. As for the case of low lying 
surface states, hot electrons in field emission resonances 
are expected to have a lifetime long enough to be scat- 
tered by steps and defects at the surface, giving rise to a 
characteristic standing wave patterns. In Fig. la a con- 
stant current dl/dV map shows clear 2-D wave patterns 
around steps and point defects on the Cu(100) surface 
(Fig. lb). The image is measured with an applied bias 
voltage (V s of 4.9 V, corresponding to the position of 
the 1st FER peak (Fig. lc). A monotonic change to 
shorter wavelengths with the applied bias (Figs. Id- If) 
reflects the energy dispersion of these states. At a first 
glance the dl/dV maps seem similar to those taken on the 
(111) faces of noble metals at lower bias. However, an 
important difference becomes apparent when looking at 
their two-dimensional (2D) Fourier Transformation (FT) 
(Figs, lg-li): here the electron wave vector is not con- 
strained to one single k(E) value but shows a broad dis- 
tribution, causin g th at the 2D FT maps resemble a disk 
instead of a ring [lfj . 

We exclude that a surface-projected bulk band, instead 
of a two dimensional state, is responsible of this broad 
distribution of k values in the 2D FT maps. Projected 
bulk bands' interference patterns might appear only at 
close distances to the sample and with oscillations cor- 
responding only to wave vectors at the band edge 
Instead, the broad distribution of parallel momentum has 
its origin in the spatial variation of the electric field along 
the surface due to the finite curvature of the tip. The lo- 
cal shift of the surface potential induced by the STM 
tip vanishes gradually with the distance away from the 
tip position (Fig. 2a). Accordingly, for a given electron 



energy E, the kinetic component along the surface di- 
rections (Ek=E-F,FER, where Efer is the FER binding- 
energy) increases continuously as the electron is accel- 
erated away from the tip. Interference patterns carry 
information of such inhomogeneous potential by show- 
ing oscillations with shorter wavelength as the tip moves 
away from the scattering potential, and a FT map with 
non-constant wave-vector distribution. Deviations from 
a perfect spherical shape of the tip apex structure are 
probably responsible for the elliptical shape of the con- 
tours. 

The evolution of dl/dV maps of Figs, ld-lf and their 
FT images show some interesting behavior. First, Fig. Id 
shows clear dl/dV oscillations with relatively short wave- 
length although the image is taken with eV s < Efer, 
i.e. below the FER onset (as shall be shown later in 
subsection , the peak position in constant current dl/dV 
spectra like in Fig. lc fits closely with the energy Efer)- 
Second, for eV s values above Efer the FT maps show 
an internal circle with a radius k m i„ increasing with the 
energy (Figs, lh and li). The origin of these phenomena 
can be understood by plotting the electron states of a 
one dimensional quantum box with a decaying potential 
as in Fig. 2a. In Figs. 2b we consider the case of an 
electronic state with energy lower than Efer- In this 
case, no allowed state exist locally under the STM tip. 
Only at distances beyond the point where the energy 
E=Efer. traveling states can exist. In the STM data 
(Fig. Id) this translates into large dark circles at scat- 
tering points, and the onset of standing waves patterns 
beyond a certain distance from the defects. In recipro- 
cal space (Fig. 2e) a continuous window of wave vector 
values from k||(0) = to kn(s) reflects the continuous 
acceleration of the electron waves from the tip position 
r|| = to the scattering point r| — s at the surface. Fol- 
lowing the same ID model, we expect that for electron 
energy above the resonance onset Efer (Figs. 2c and 2d) 
a minimum wave vector k||(0) = 2m*/H 2 (E — Efer) 1 ^ 2 
appears in the ID plots (Figs. 2f and 2g), corresponding 
to the internal circle in Figs, lh and li. 

Therefore, the width of the wave vector distribution 
will reflect the spatial change of the local (tip-induced) 
potential shift. Ideally, for a given value of the tip radius, 
the critical angle determined by the exponential decay of 
the Fowler-Nordheim transmittivity 27|, permits to es- 
timate the maximum parallel component of the electric 
field at the tip and, therefore, the change of parallel mo- 
mentum, based on simple classical trajectory considera- 
tions. Assuming a radius of curvature of the tip R=10 nm 
and a tip-sample distance Z=15 A we estimate a change 
in parallel momentum of 2 nm _1 at 5 V, in agreement 
with our previous analysis shown in figures lg)-i) . Re- 
sults similar to those shown in Fig. 1 are observed for 
other FERs that appear below 8 eV. However, no wave 
patterns are seen above this value, indicating a significant 
decrease of the electron lifetime. 



FIG. 1: (Color online) a) dl/dV map of the Cu(100) surface shown in b) (1=10 nA; Vs=4.9 V). Standing waves are clearly 
seen around steps, surface and sub-surface defects, c) Constant current dl/dV spectra around the first field emission resonance, 
d-f) dl/dV maps of a smaller region taken at the bias values indicated in the figure, g-i) 2-D Fourier transform of d-f). The 
arrows indicate the inner and outside radius of the disk appearing in K-space (see the text). 



Line shape of dl/dV spectra 

To understand the role of the surface electronic struc- 
ture in the dynamics of FERs, next we explore the in- 
formation contained in dl/dV spectra about the energy 
width of FERs by analyzing their line shape and com- 
pare with lifetime estimations based on a phase coherence 
length analysis of interference wave patterns [28| . 

In Fig. 3 we show a series of constant current dl/dV 
spectra at different set point current values. Sharp peaks 
appear in a wide sample bias range covering the en- 
ergy range of the projected Cu bulk gap and above [29| . 
Each peak corresponds to the onset of a 2-D resonance 
state. In field emission regime, the tunnel transmit- 
tivity T(E,V) is sharply peaked at the tip Fermi level 
(E = Ef); most of the current comes from a narrow 
energy window of about 100 meV below Ef (22], caus- 
ing FERs to appear as peaks in dl/dV spectra. On the 
contrary, in case of surface states close to Ep dl/dV spec- 



tra show a line-shape close to a step 14]. Interestingly, 
we find that the apparent width of the dl/dV peaks ex- 
hibits a non-monotonic behavior (shown in the inset). 
Resonances lying inside the gap are narrower than those 
appearing above approximately 8 eV, whose width in- 
creases in agreement with the steep decrease of surface 
reflectivity as the top edge of the projected bulk band 
gap is crossed [1, [29| . 

The finite line-width of resonance states in the gap is 
expected to be dominated by intrinsic factors like the 
electron reflectivity of the surface, spatial extension of 
the wave function 13 ll and electric field strength at the 
tip-sample region [301 ] - However, as it is shown in the 
inset of Fig. 3, the apparent width of the first three 
peaks in constant current dl/dV spectra is considerably 
larger than the intrinsic line-width estimates based on a 
phase coherence length analysis [28j|, which for the first 
FER agree with previous ab initio calculations [30l | and 
estimations .15] 

for the same system. Therefore, an ad- 
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FIG. 2: (Color online) a) Schematic representation of the non-conservation of parallel momentum due to the applied field by a 
tip of finite radius R. At a given value of electron energy E — eV s — 4> the kinetic energy for parallel motion £V. increases from 
the tip position rn = to distant sites at the surface r| — s. b), c) and d) show the charge density of electron states confined 
in a large box with a potential V(x) — Vb/[1 + (x/xo) 2 ] with energy below, slightly above and above Vb, respectively [Vb = 0.03 
(arb. units) and xo = 200 (arb. units)], e), f) and g) show the corresponding FT. Curves c), d), f) and g) are shifted vertically 
for clarity. The arrows indicate the range of wave vectors with significant weight in the FT and, thus, confirm the measured 
behavior in Figs. Ig), fh) and fi). 



ditional broadening mechanism must exist to explain the 
apparent width values in these dl/dV plots, which, pre- 
sumably, is related with the method of measurement. 

During the acquisition of dynamic (constant current) 
dl/dV spectra, the distance (Z) vs. bias voltage (V) char- 
acteristics (Z(V)) exhibit a pronounced step- like behavior 
as the resonance is crossed for electron energies lying in 
the projected gap. The resonant electron transmitivity in 
this close feed-back loop spectroscopy is expected to be 
affected by the continuous change of the tunneling bar- 
rier shape with bias voltage and tip-sample distance. It is 
reasonable to assume that this dynamic mode will intro- 
duce some distortion in the resonance's line shape respect 
the ideal static situation, in which the tunneling barrier 
shape is kept fixed at Z(V) peak values. It is only in this 
latter case that one could relate the width of peaks in the 
transmitivity T(E) to the intrinsic width of resonances. 
We have performed a model calculation to establish a link 
between intrinsic energy line-width of resonances in the 
static tunneling transmitivity T(E) and the correspond- 
ing apparent width in constant current dl/dV spectra, 



which is related to the auxiliary dynamic transmitivity 
T(E F ,V) at tip the Fermi level. 

Our calculations are based on 1-D model potentials for 
the tip and sample including the work function and Fermi 
energy that defines the bottom of the surface potential. 
For the surface, a periodic sinusoidal modulation that de- 
termines the magnitude and position of the energy gap 
at the r point of the surface Brillouin zone 32] is in- 
cluded. Parallel wavevector components are considered 
in the full 3-D calculations assuming a free-electron-like 
(parabolic) dispersion. To model electron transmission 
inside the gap, inelastic scattering at the surface is in- 
cluded by usin g a complex potential, similar to previous 



LEED studies [33j, |3J]. Its imaginary part (Vbi) intro- 
duces a decay of the electron flux due to absorption. A 
smooth matching between the tip and sample potentials, 
which includes multiple images, is used 35]. The cal- 
culation of the tunneling current for a given tip-sample 
distance (Z) and bias voltage requires the knowledge of 
the energy dependence of the barrier transmission co- 
efficient T(E,V) below the tip's Fermi level Ep (3(|. A 
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Sample bias (V) 

FIG. 3: (color online) dl/dV spectra taken for several current set points ranging from 0.6 nA to 60 nA (V ac =10 mV rms). The 
inset shows the width of the peaks from a Lorentzian fit, as well as values obtained from a phase coherence length analysis [2^ ]. 
Shaded grey areas correspond to the energy location in F of copper bulk bands projected on the (100) surface [291 ]. 
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FIG. 4: (color online) a) Comparison between a calculated dl/dV spectrum (black curve) and the transmission coefficient (red 
curve) at the tip Fermi level T(_B_f;V). The inset shows a comparison between a static calculation (black) of T(E) at Z(V) values 
of the first resonance (thus reflecting the intrinsic FER line shape), T(i5_F;V) (red) and dI/dV(V) (blue). b)dI/dV spectra for 
two Voi values (see the text). 
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quantitative agreement between measured and calculated 
dl/dV spectra is not persecuted in this model approach, 
mostly due to the lack of knowledge of effective tunneling 
areas, tip-sample distance or tip work function. Instead, 
we can provide a qualitative picture of the effect of the 
dynamic measurements on the peaks' width. The con- 
stant current dl/dV spectra are calculated by numerical 
differentiation of the current I(V,Z) along the constant 
current Z(V) characteristic. 

In Fig. 4a we show a comparison between calcu- 
lated dI/dV(V) spectra and the corresponding dynamic 
transmitivity T(£t?;V). Both curves exhibit a similar 
shape and a characteristic increase of their line-width 
with applied bias. This confirms the high collima- 
tion of field emitted electrons in a narrow energy win- 
dow below Ep. The inset compares the conductance 
dI/dV(V) around the first peak, the dynamic transmi- 
tivity T(Ep-,V) and the intrinsic line-shape of the corre- 
sponding FER obtained in a static calculation at Z(V) 
peak values. Lorentzian fits to T(Ei?,V) and dI/dV(V) 
give width values of ~150 meV and ~350 meV, respec- 
tively, while the intrinsic width of the resonance in T(E) 
is ^100 meV. Therefore, the broader line shape in dl/dV 
spectra must be related both to the above mentioned fi- 
nite energy collimation and the variation of the tunneling 
barrier shape with applied bias. For FERs in the gap, the 
increase in electric field strength with bias shifts the res- 
onances to higher energy, appearing as broader peaks in 
dI/dV(V). The increase of tip-sample distance in the dy- 
namic method of measurement partially compensates for 
this broadening effect. The best conditions for a quan- 
titative line-shape analysis can be achieved at constant 
field strength conditions (~V/Z) and low set point cur- 
rent values (low applied field). 

By comparing (Fig. 4b) the shape of dl/dV curves cal- 
culated for two different values of Vbi we find that only 
the first peak in the gap broadens as a response to the in- 
crease in absorption (inelastic scattering). This confirms 
that inelastic scattering is the broadening mechanisms 
of FERs lying inside the projected band gap. At ener- 
gies above the gap inelastic effects play a minor role 37 
and the FER intrinsic line-width is dominated by elastic 
coupling to the bulk continuum (shaded area). In this re- 
gion, the intrinsic resonances' line-width is considerably 
larger (hence, no wave patterns in dl/dV maps could be 
seen for this energy) and, also, V/Z is almost constant. 
It is then expected that the effect of the dynamic broad- 
ening will be smaller, and the experimental dl/dV peaks' 
apparent line-width will be close to the intrinsic value. 



CONCLUSIONS 

We find that the scattering of electrons in field emis- 
sion states with defects in a Cu(100) surface gives rise to 
isotropic standing wave patterns, which reflect their dy- 



namics in response to the electric field gradients induced 
by the STM tip at the tip-sample interface. Through 
a combined theoretical and experimental study we have 
identified: (i) the origin of characteristic peaks width in 
constant current dl/dV spectra above the vacuum level as 
a combination of both, their FER intrinsic line shape and 
the extrinsic distortion due to the measuring process, and 
(ii) for FERs in the gap this distortion introduces an ad- 
ditional broadening leading to a non-monotonic behavior 
of the width with sample bias. Our results show that STS 
in the field emission regime can be used to gain informa- 
tion about the electron dynamics and surface electronic 
properties at energies well above the vacuum level. 
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